Antiferromagnetic correlations and impurity broadening of NMR linewidths in 

cuprate superconductors 
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We study a model of a d-wave superconductor with strong potential scatterers in the presence 
of antiferromagnetic correlations and apply it to experimental nuclear magnetic resonance (NMR) 
results on Zn impurities in the superconducting state of YBa2Cu307_5 (YBCO). We then focus on 
the contribution of impurity-induced paramagnetic moments, with Hubbard correlations in the host 
system accounted for in Hartree approximation. We show that local magnetism around individual 
impurities broadens the line, but quasiparticle interference between impurity states plays an impor- 
tant role in smearing out impurity satellite peaks. The model, together with estimates of vortex 
lattice effects, provides a semi-quantitative description of the impurity concentration dependence of 
the NMR line shape in the superconducting state, and gives a qualitative description of the tem- 
perature dependence of the line asymmetry. We argue that impurity-induced paramagnetism and 
resonant local density of states effects are both necessary to explain existing experiments. 

PACS numbers: 74.25.Ha,74.72.Bk,74.81.-g 



I. INTRODUCTION 

The response of a correlated electron system to a lo- 
cal perturbation can often provide important information 
about the ground state of the pure system. This princi- 
ple has been successfully applied to the substitution of 
various impurities, particularly Zn, Ni, and Li for Cu 
in the Cu02 planes of various high-Tc superconductors, 
particularly in YBa2Cu307_5 . In addition to traditional 
studies of the effect of impurities on bulk properties, local 
probes like NMR and scanning tunnelling spectroscopy 
(STS) have provided considerable information about how 
electronic wavefunctions are distorted near the impu- 
rity site. In the normal state of this system, NMR has 
shown that irnpurities enhance local antiferromagnetic 
correlations P, and in the presence of the applied DC 
field display a staggered pattern of magnetization which 
decays over a few lattice spacings. This polarizability Sx 
has moreover a characteristic temperature dependence, 
which is Curie-like Sx in the underdoped system, 

but evolves to Curie- Weiss-like behavior 5x ~ [T + 9)~^ 
in the optimally-to-overdoped range 0. Because in- 
creases rapidly with doping, it has sometimes been inter- 
preted as a Kondo temperature, enhanced in the presence 
of higher carrier densities capable of screening the mag- 
netic moment induced by the impurity. Other pictures 
of this phenomenon, which do not rely on Kondo screen- 
ing, have been put forward as well. For example, in the 
weak-coupling approaches of Bulut|l,|3 and OhashiQ, 
an extended potential is found to produce a Curie- Weiss- 
like local susceptibility in a Hubbard model treated in 
mean field, due to the coupling of the antiferromagnetic 
q = (tt, tt) response of the lattice system to the uniform 
q = response by the inhomogeneity. 

In the superconducting state, interpretation of the 
NMR signal is complicated by the intrinsic field distri- 



bution introduced by the vortex lattice, and by the van- 
ishing of the q = susceptibility in the singlet pair state. 
On the other hand, OhashiQ argued that the mode- 
coupling effect induced by the inhomogeneity persists and 
is relatively enhanced by the opening of the gap. Experi- 
mentally, the enhancement of the local susceptibility was 
indeed found below Tc!^, and tended to a large constant 
value at very low T. Recently, Ouazi et ai. 9] measured 
the evolution of the ^^O NMR line with increasing Zn 
concentration, and observed the formation of the stag- 
gered polarization cloud for the first time in the super- 
conducting state. They argued that the primary line shift 
was due to the nearly field-independent vortex distribu- 
tion, and that the broadening was a combination of the 
enhanced impurity effect and the simultaneous narrow- 
ing of the vortex field distribution due to the increased 
penetration depths A in the dirtier systems. Missing from 
this picture is an understanding of the magnitude of the 
impurity broadening and how it really occurs; if one con- 
siders only a single impurity one expects large values of 
the magnetization on the Zn nearest neighbor sites, which 
should lead to a well-defined satellite line as in the case of 
NMR on Li impurities in the normal state 3] . These have 
not been detected in samples with Zn concentrations at 
the per cent level. 

An alternate picture of the observed phenomena in the 
superconducting state is obtained if one considers the lo- 
cal susceptibility due to quasiparticles in the d-wave su- 
perconductor, proportional to the local density of quasi- 
particle states (LDOS) at the Fermi level. A significant 
enhancement of the local susceptibility is then to be ex- 
pected from quasiparticle bound states alone. Williams 
et aJ.^3l proposed that these quasiparticle resonant 
states-corresponding to those imaged by STM experi- 
ments around Zn atoms in Bi2Sr2CaCu208+2; (BSCCO- 
2212)[ll]-might be entirely responsible for the enhanced 
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magnetic response near Zn seen in NMR. Chang et aJ.|l2| 
then argued that for a single nonmagnetic impurity, the 
temperature dependence of the observed spin-lattice re- 
laxation time and Knight shift could indeed be qualita- 
tively understood in terms of LDOS enhancement due to 
impurity bound states alone. 



There are some difficulties with the naive interpreta- 
tion of the NMR measurements entirely in terms of the 
LDOS enhancement near impurities, however. First of 
all, significant T-dependent enhancements of local sus- 
ceptibilities near nonmagnetic impurities occur in the 
normal state of optimally doped cuprates as well. An 
"LD OS-only" approach cannot account for this since im- 
purities do not produce LDOS resonances in the nor- 
mal (metallic) state. Secondly, the NMR experiments 
on these materials clearly show that the magnetization 
near a Zn alternates in sign. This is incompatible with 
the paramagnetic character of the quasiparticle Pauli sus- 
ceptibility, i.e. the susceptibility enhancement due to im- 
purity bound states is always positive, so while Friedel- 
type oscillations can occur, the magnetization is always 
aligned with the external field. Finally, the "Knight 
shift" calculated in Ref. is defined to be a local suscep- 
tibility enhancement very near the impurity; in fact the 
measured Knight shift in experiments where the nucleus 
is distinct from the impurity itself is the shift of the total 
NMR line, determined by sites far from the impurities. 



Thus a theoretical calculation which includes local 
magnetic moment formation, together with quasiparticle 
impurity bound states and their interference, is of con- 
siderable interest to understand the simple but striking 
features of the NMR experiments in the superconduct- 
ing state 8, 9J. A complete theory of these phenomena 
must be able to account not only for the Knight shift, 
but for the detailed behavior of the site-specific NMR 
lines produced by the different nuclei probed in different 
experiments. In this paper, we study Zn ions modelled 
as strong potential scatterers in a d-wave superconduc- 
tor with AF correlations treated within a weak coupling 
approach as in Ref. 0- In Section ^ we describe the 
model for a single impurity in a d-wave superconductor 
with correlations and the magnetization it induces, then 
in Section IIIII study interference effects on the magne- 
tization distribution when many impurities are present. 
In Section [l VI we combine the predicted impurity contri- 
bution to the NMR linewidths with vortex effects, and 
compare to the Ouazi et "'^'^O NMR experiment. Sec- 
tion is devoted to the application of the same results 
to compare with experiments using the ^Li nucleus. In 
Section IVII we present our conclusions and implications 
for other experiments on the cuprates, as well as new 
questions raised by our interpretation. 



II. SINGLE IMPURITY IN SYSTEM WITH 
ANTIFERROMAGNETIC CORRELATIONS 

A. Formalism 

A strong nonmagnetic impurity in the presence of 
antiferromagnetic correlations will induce a pattern of 
local staggered magnetization with maximum peaks 
on nearest-neighbor sites. When many impurities are 
present, these local states interfere, washing out nearest- 
neighbor peaks and producing a smooth distribution of 
local fields. The NMR line shift and broadening is caused 
by this impurity effect in conjunction with the vortex con- 
tribution. In order to study this magnetic behavior, we 
begin with a two-dimensional tight-binding Hamiltonian 
of a d-wave superconductor with AF correlations treated 
within mean field theory: 

ij(T i(7 

+ E(^^44^i+H.c.), (1) 

iS 

where the hopping term includes nearest-neighbor hop- 
ping t and next-nearest-neighbor hopping t' , eig. = Vimp ~ 
gfXB^Ba describes the impurity and Zeeman site en- 
ergies, ^ is the chemical potential, A^.; is the near- 
est neighbor pairing potential, and 5 G {x, — x, — y} 
are unit lattice vectors to nearest- neighbors. Here 
g K, 2 \s the electron g- factor, jiB is the Bohr mag- 
neton, and B is the applied field along the c-axis. 
The electron number and d-wave (singlet) pairing pa- 
rameters are defined as rijo- — {c\^Cic) and A^^ — 
V {ci^Ci+si + Ci+5f Ci|) /2. We note that the Hamiltonian 
has been used extensively to study bulk competing 
phases, disorder and vortex induced magnetization, as 
well as novel bound states at interfaces between antifer- 
romagnets and superconductors|l3j . 

Eq. (P) can be diagonalized by using the Bogoli- 
ubov transformation. The corresponding Bogoliubov-de 
Gennes equations must be solved iteratively until a self- 
consistent solution is found: 

where positive eigenvalues correspond to spin-up ex- 
citations and negative eigenvalues correspond to spin- 
down excitations. The matrix operators are defined 
by icUns = - Y.ij iij'^nj + {Uui^^ + t^cr - [i.) u„,i and 
Aun^i = AsiUn.i+s- The mean field parameters, up- 
dated after each iteration until sufficient convergence is 
achieved, can be computed by: 

n \ B / 
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and riia = {rii + ami) /2, where rii is the average electron 
density at site i and is the magnetization on site i: 
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It is important to note that in this model the supercon- 
ducting pairing is taken to be a phenomenological con- 
stant, and is not affected by the Hubbard repulsion U. 
In this sense results may differ from true inhomogencous 
spin fluctuation models where both magnetism and pair- 
ing are driven by the same correlations. 

When deciding which ordered phase is the ground 
state, it is necessary to know the total energy E of the 
system, defined as (H). It is given by: 



E = K.E. + 
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where K.E. is the kinetic energy of the system 
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K.E. = - u«j<,j] tanh (^^j -(7) 
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FIG. 1: T = U-n phase diagram for Hubbard plus nearest- 
neighbor pairing model treated in mean field used in this 
work. F, AF and P are the usual ferromagnetic, antiferro- 
magnetic, and paramagnetic phases found for the homoge- 
neous Hubbard model. LAF is the local antiferromagnetic 
phase where spontaneous staggered moments form around a 
single potential scatterer. The gray lines indicate the bound- 
aries of coexistence of d- and s-wave superconductivity with 
the magnetic phases. Calculations are performed in the dark 
rectangle region between < U < 2.0 and at 15% hole doping, 
n = 0.85. 



We will often speak of the total magnetic moment (or 
spin) of the system. This is given (in units of h) by 

S.^^Y."^,. (8) 

i 

We will also make reference to the d-wave order param- 
eter, defined as 

d^ = ^ [A^,, + A-^,^ - Ay^, - A_y^,] . (9) 

In the following we operate at a band filling n = 0.85 
such that the ground state of the homogeneous system is 
always pure d-wave of the form At oc cos — cos ky . 

Finally, we have occasional need for spatially resolved 
spectral information, and so also calculate the LDOS 
N{E,i), via 

NiE,i)= i^(|u„,/<5(£;-^„) + |i;„,.|'<5(^ + ^„)) . 

n 

(10) 

The chemical potential fi is adjusted to produce an aver- 
age electron density of 0.85, corresponding to 15% hole 
doping characteristic of optimally doped cuprates. We 
give all energies in units of t and set t' = —0.2, to mimic 
typical Fermi surface shapes found in these systems, and 
y = 1 to give a critical temperature of 0.15t. 



B. Single impurity phase diagram 

In Fig. ^ we give a schematic phase diagram for the 
model with crude phase boundary lines and details 
suppressed for simplicity. Here, "AF" , 'F" , and "P" de- 
note self-consistent mean-field phases of Eq. char- 
acterized roughly as antiferromagnetic, ferromagnetic, 
and paramagnetic, respectively. For example within the 
"AF" region it is well-known that spin density waves 
with ordering vectors other than exactly Q = (tt, tt) can 
be stabilized with details depending on the doping level 
and band-structure parameters. Since we have included 
a separate nearest-neighbor pairing interaction term V 
in the Hamiltonian in order to study the superconduct- 
ing state, we have also indicated in Fig. ^ the regions 
of doping over which nearest neighbor d-wave or s-wave 
pairing symmetry characterizes the ground state. Note 
that in this case the ordered magnetic phases coexist with 
superconductivity within the model; again details have 
been suppressed for simplicity, and because we are here 
primarily concerned with the paramagnetic (d-wave su- 
perconducting) phase. 

Upon addition of a single strong impurity potential 
to the model, we find new inhomogeneous ground states 
present. Of most interest is a region of local staggered 
magnetism surrounding the impurity, referred to as the 
"local antiferromagnetic" (LAF) phase. In this phase. 
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FIG. 2: Top: phase diagram for a single impurity as a function 
of the impurity potential Vimp and U, as determined by the 
presence of a nonzero magnetization at T = 0.013. Note 
there are only two phases, ^^^O and Sz ~ 1/2. Bottom: 
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100, dependence on U of 



the staggered magnetization, defined as Mi = (— l)*mi, for 
a 28 X 28 system. The system is completely nonmagnetic 
{M = 0) for small values of U, assumes a local staggered 
state for intermediate values, and saturates toward a roughly 
homogeneous AF phase for large values. For the cases with 
nonzero magnetization, Sz = 1/2. The wave-like AF ordering 
in the bulk for U = 2.5 and 4.0 is due to finite size effects and 
the periodic boundary conditions of the system. 




FIG. 3: d-wave order parameter suppression at an impurity 
of strength Vimp = 100 for [/ = 1.75, V = 1, 34x34 system, 
giiBB/2= 0.004, and T =0.013. 
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the impurity-induced staggered magnetization vanishes 
at a large distance from the impurity, and the net spin 
Sz summed over the whole system is found to be l/2|3l- 
At larger [/, the impurity still generates a net spin 1/2, 
but the long-range ordering dominates, and the staggered 
magnetization has the maximum polarization arbitrarily 
far from the impurity. For a given fixed band struc- 
ture and doping n — 0.85, we have plotted in Fig. |21 
the transition line between the (total) Sz = state and 
the Sz — 1/2 state for varying U and impurity potential 
Vimp] clearly increasing either U or Vimp tends to favor 
the local magnetic "spontaneous moment" state. In the 
lower panel of Fig. [21 we show the staggered magnetiza- 
tion patterns in each of the states for increasing U. Note 
that in general we do not find an intermediate state with 
S — and impurity-induced magnetism, in contrast to 
a recent study of the same Hamiltonian in Ref. Q, and 
believe the presence of this state to be an artifact of the 
particular size system studied by these authors. We find 
generically either a S'z = 1/2 state with local staggered 
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FIG. 4: The LDOS for a single-impurity system dis- 
playing paramagnetic behavior {U = 1.25, T = 0.013, 
gfiBB/2=0.004), computed using a 28x 28 system and 20 
X 20 supercell. This result is in agreement with past studies 
using a similar model 



magnetic order, or a state with no local magnetization 
at all in zero field, referred to a,s Sz = in Fig. |21 This 
is consistent with the results of Wang and Lee 15j, albeit 
derived in the weak coupling limit rather than for the 
t — J model. 

We recover other single-impurity results known for this 
model which we state for the sake of completeness. For 
example, the d-wave order parameter is strongly sup- 
pressed in the vicinity of a strong potential scatterer, 
mostly over a length scale of one lattice spacing, but with 
a longer-range decay envelope over the coherence length 
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FIG. 5: The field-induced local magnetization for U — 1.75, 
giiBB/2 = 0.004 and T = 0.013 on a 34x34 system. The total 
moment of the system is {S^) — 0.294. Weak long-range AF 
correlations extend out from the impurity along 45° diagonals. 



^0 (Fig. EJ- The impurity also induces a much smaller 
s-symmetry order parameter(not shown), and gives rise 
to a strong resonance in the LDOS near the Fermi level. 
The spatial intensity of the resonant state is centered pri- 
marily on the nearest-neighbor sites of the impurity (Fig. 
01, in apparent conflict with the most naive interpreta- 
tion of STM experiments; there are several competing 
explanations why this could be so. These phenomena 
have been reviewed and references given in jig (see also 



C. Response of Sz = state to applied field 

In this paper we focus primarily on the impurity- 
induced state which has net spin = in zero field. 
This is because there is no evidence for impurity-induced 
magnetization of any kind in zero field for almost all of 
the YBCO phase diagram J^. NMR experiments, which 
are sometimes cited as providing evidence for sponta- 
neous impurity induced magnetization, are of necessity 
performed in finite applied field. Current B = neu- 
tron scattering 19] and /iSR|20|| measurements find no 
evidence of ordered static magnetization at any wavevec- 
tor. Furthermore, both NMR and direct susceptibility 
measurements indicate that the induced states are para- 
magnetic, i.e. the magnetization vanishes proportional 
to the applied field. We note that there is considerable 
recent evidence that the situation is different in LSCO, 
where static magnetism appears to exist even without an 
applied external field at low temperatures [2]1 12^. This 
may also explain unusual transport properties in LSCO 
compared to YBCO 23, 24]. 

With these considerations in mind, we choose a value 
of U which will induce significant antiferromagnetic cor- 



relations close to half-filling, but is not sufficient to cause 
the formation of magnetic moments around impurities in 
zero field. The impurity is taken to have an on-site spin- 
independent potential of strength Vimp = 100, roughly 
consistent with STM at least as far as the energy of the 
Zn LDOS resonance is concerned. To start, we study the 
Bogoliubov-de Gennes equations for a single such impu- 
rity with systems of size 34 x 34 with periodic bound- 
ary conditions. The Zeeman response of the electronic 
spins to the applied field is included in the eta term in 
Eq. (P), as discussed above. In our approach, we do 
not include the orbital response of the system to the ap- 
plied vector potential, but rather attempt to account for 
the presence of the vortex state phenomenologically (see 
below). We wish to compare our results with the data 
from the Ouazi et experiment; we therefore take 

t = 100 meV and set \gnBB/2\ = 0.004 (5 « 7 T) 
and ksT = 0.013 (T 15 K ). Application of such a 
field induces, as expectedP, 0, |E Hi a local staggered 
magnetization of the Cu spins around the impurity site, 
depicted in Fig. \S\ 

Below we argue that results consistent with experiment 
on optimally doped samples require a value of U close to 
the threshold for the creation of static impurity-induced 
zero-field magnetism (i.e. close to the phase boundary 
in Fig. [2J|. This means that the field dependence can 
acquire nonlinearities, as shown in Fig. For most of 
this work, we examine U = 1.75, although we also exhibit 
the consequences of choosing other values. 

The size of t assumed in order to compare with exper- 
iment is probably a factor of 2-4 smaller than deduced 
from ARPES experiments on cuprates, but we do not 
expect this to alter our qualitative conclusions. The low 
value of t is chosen such that reasonably low tempera- 
tures T/Tc can be accessed without encountering finite 
size effects. For this choice, together with the choice of 
pair interaction V, Tc = 0.15 corresponds to about 175K, 
and the field parameter of gfj,B/2 = 0.004 to the exper- 
imental value of 7 T. The reader should therefore take 
the scales given in Kelvin only as extremely rough com- 
parisons. 

In the linear field response regime (Fig. 1^, our calcu- 
lations should be very similar to those of Ohashi for the 
nearest-neighbor weak-field susceptibility|25. Thus it is 
not surprising that we also find a strong increase of the 
magnetization on the nearest neighbor sites as the tem- 
perature is lowered. This increase is weak in the normal 
state, then slows slightly at the superconducting tran- 
sition as the gap opens, as depicted in Fig. [7| As the 
temperature is lowered further, the resonant state in the 
d-wave superconductor forms (Fig. driving the sus- 
ceptibility to a large value, which however is expected to 
saturate at T ^ 0, as indicated in the Figure. We are un- 
able to calculate results accurately below a temperature 
T ~ D{a/L)'^, where D w 8i is the bandwidth, because 
the thermal energy becomes of order the level spacing 
in the finite size system. For the simulations reported 
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FIG. 6: The response Sz induced by a field B of a system 
with U — 1.75 for three values of T. In the limit of zero field, 
the system is nonmagnetic regardless of temperature. Inset: 
field response for U — 1.25. 
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FIG. 8: Top: magnetization along (1,0) direction through 
impurity with on-site potential Vimp = 100, g^BB/2 — 0.004, 
U = 0,1, 1.5, T = 0.013. Bottom: same but along the (1,1) 
direction. 
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FIG. 7: Magnetization in fixed field g/iflB/2=0.004 vs. T on 
nearest neighbor and next nearest neighbor sites for a strong 
impurity with U = 1.75. 



est neighbor susceptibility is also enhanced, but has the 
opposite sign because the correlations are antiferromag- 
netic. 

It is worth noting that, in the presence of the nonzero 
external field, a modulated local magnetic state is present 
even in the absence of the antiferromagnetic correlations 
driven by U. This is the analog, in the d-wave super- 
conducting state, of Friedel-like spin density oscillations 
which represent the response of the normal metal to a 
local perturbation. As such, the oscillations necessarily 
take place at an incommensurate wave vector 2fci? which 
is however close to (tt/o, tt/o) because the system is close 



to half-filling. Note that these U — magnetization oscil- 
lations, shown in Fig. |S1 are driven by the Pauli suscepti- 
bility in the d-wave superconductor. For a weak impurity, 
this response is quite weak at low temperatures, due to 
the linear w— dependence of the d-wave density of states 
near the Fermi level. On the other hand, the LDOS reso- 
nance at the Fermi level in the case of a strong impurity 
enhances this local response substantially. Note that the 
magnetization is always positive, however, since the lo- 
cal susceptibility is proportional to an enhanced LDOS 
at the Fermi level. When correlations are added, as in- 
dicated by the increasing U in Fig. |S1 the response can 
be many times that of the pure BCS system with non- 
interacting quasiparticles, and takes on an alternating 
character, as seen. These effects will lead to asymme- 
tries in NMR lineshapes, as discussed below. We note 
further that the "background" homogeneous magnetiza- 
tion of the system in nonzero external field is present in 
Fig. |S1 but barely visible due to the small value of the 
homogeneous d-wave susceptibility at low T. 



7 



c 
o 



N 




FIG. 9: Magnetization m for one configuration of a system 
with 1.5% random impurities and U = 1.25, T = 0.013, B = 
0.004. 



III. MANY-IMPURITY MAGNETIZATION 

In the presence of many strong impurities, the wave 
functions of electrons bound to the impurity interfere at 
long distances, leading to collective behavior which is no 
longer describable by the 1-impurity model. These effects 
have been studied in d-wave superconductors without an- 
tiferromagnetic correlations HI El HI S El El 
The interference of the many impurity states leads to a 
splitting of bound state energies and an accumulation 
of low-energy impurity-induced energy eigenvalues which 
are spread out over a so-called "impurity band" . In the 
(i-wave case, the formation of the impurity band and 
the corresponding quasiparticle localization problem are 
strongly influenced by the fact that significant overlaps 
between two impurity states can take place only if the im- 
purities are "oriented" with respect to one another such 
that the nodal quasiparticle wavefunctions overlap along 
the 110 direction. In fact, analysis of the two-impurity 
problem show that interference effects can take place over 
many tens of lattice spacings between optimally oriented 
impurities I2E El E3|- systems with per cent level dis- 
order, however, these effects are significant also for pairs 
of impurities aligned along the 100 direction. 

In the presence of correlations, interference effects ex- 
tend to the magnetic channel and are enhanced by in- 
creasing U. In Fig. 13 we show a system in applied field 
in the presence of many strong impurities. It is clear 
that the size of the magnetization on nearest neighbor 
sites varies significantly according to the local disorder 
environment. To clarify this, we compare a few of these 
impurities in the many-impurity sample with the compa- 
rable impurity in isolation in Fig. [TUl 

The distribution of magnetizations shown in Fig. ^ 
represents all the information necessary to calculate the 
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FIG. 10: 0° magnetization cuts through 1% impurity sys- 
tems (a-c)and a single isolated impurity (d), taking U = 1.25, 
T = 0.013, gyLBB/2 = 0.004. Arrows indicate a cut passes 
through an impurity site, where magnetization is zero, and 
large peaks occur on nearest-neighbor sites. Interference ef- 
fects can both enhance and suppress nearest-neighbor mag- 
netizations relative to the 1-impurity result shown at bottom. 
Note the cut in panel (c) does not pass through an impurity 
site, but passes through two nearest-neighbor sites. 



NMR response within our approach, since it determines 
the distribution of nuclear spin precession frequencies, in 
different linear combinations depending on the locations 
of the nuclei relative to the Cu sites. We focus in more de- 
tail on the actual NMR lines in Sections ll VI and Ivl below . 
For now, we are interested in showing how interference ef- 
fects influence the "bare" distribution of magnetizations. 
To this end, we collect all the magnetizations in the sys- 
tem in a histogram for a single impurity, then compare 
to progressively larger impurity concentrations in Fig. 
im For a single impurity, small satellite peaks are visible 
in the spectrum since the same magnetization value ap- 
pears on all sites with fourfold symmetry. Note that the 
satellite peaks associated with the nearest neighbor sites 
occur at magnetization values well outside the range of 
the plot! With the addition of a few random impurities, 
these magnetization values are split, the distribution is 
smeared, and satellites are seen to disappear already at 
sub-percent level concentrations. In addition, positive 
magnetizations are seen to be preferentially enhanced, 
due to the density of states effects discussed above. 
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FIG. 11: Magnetization histogram for a 34 x 34 system 
with (7=1.75, and 1,4,8 and 12 impurities at T = 0.013 and 
giJ.BB/2 = 0.004. 
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FIG. 12: (Color online) Distribution in real space of mag- 
netization values around a single Vimp = 100 impurity with 
(7=1.75, r = 0.013, and gfJ.BB/2 = 0.004. 



It is interesting to ask why interference effects in the 
magnetic channel are so important in the superconduct- 
ing state that satellite features are immediately elimi- 
nated. To this end we plot the distribution of magne- 
tization values in real space in Fig. El It is seen that 
those values which contribute to the satellites close to the 
peak (which eventually determine the width of the line 
for a finite density of scatterers) are primarily located 
in the 45° tails of the quasiparticle wavefunctions some 
10-15 lattice spacings from the impurity site (see below, 
however) . The orbital parts of these wave functions are 
known to interfere strongly in the absence of correlations 
provided other impurities are appropriately oriented, so 
it is no great surprise that the magnetic parts of these 
wavefunctions also strongly interfere. 



Pure, Normal State 
Pure, SC 
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FIG. 13: Magnetic susceptibility per site {1/L )dSz/dB for 
a 28 X 28 system at i7 = 1.25t for normal state (d = 0), 
superconducting state, one impurity and 1% impurities. 



Any bulk measurement of magnetization results in an 
average over this smeared magnetization distribution. In 
addition, the temperature dependence of the magneti- 
zation depend on their position relative to the impu- 
rity. Contributions from sites far from impurities de- 
crease with decreasing temperature, as for the homoge- 
neous d-wave superconductor. Impurity nearest-neighbor 
susceptibilities are strongly enhanced, on the other hand. 
These effects combine to determine the total temperature 
dependence of thermodynamic properties. For example, 
if one measures the total susceptibility of the sample, it 
exhibits an upturn at low T if the density of impurities 
is a significant fraction of the sample (Fig. I13|l . 



IV. IMPURITY CONTRIBUTION TO "O LINE 

A. Disorder dependence of ^^O line 

The planar ^^O nucleus is situated halfway between 
two Cu sites, or between a Cu and a Zn. It is assumed 
that it senses the local field proportional to the sum of 
the magnetizations on the two sites closest to it. For ex- 
ample, O nuclei far from the impurities are subjected to 
'^XhomB, where Xhom is the susceptibility of the homoge- 
neous system, whereas an O nucleus next to a Zn atom 
is subjected to (xnn + 0)i?, since there is effectively zero 
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electron density on the Zn site by assumption. Nuclei at 
varying distance from the impurity will measure different 
combinations of local magnetizations. 

In an NMR experiment, the resonance frequency v of 
a nucleus is given by: 

I. = -^S [1 + Xo,b + if,p,„] , (11) 
ZTT 

where 7 is the gyromagnetic ratio of the nucleus, B is 
the applied magnetic field, Ko^^, is the T-independent 
orbital contribution of valence and inner-shell electrons, 
and Kgpin is the spin contribution from the electrons. 
In a simple metal or in a cuprate above Tc, Kgpin is 
proportional to the uniform electronic spin susceptibility 
Xspin through: 

^ spin — 5 V J 

where Ah / is the hyperfine coupling between the nucleus 
and the electrons, and fis the Bohr magneton. In the 
specific case of ^^O NMR for planar oxygens, a ^^O 
nucleus at position (x; y) is coupled to the spin suscep- 
tibility through its two neighboring Cu, leading tojS^I : 




-100 100 200 300 



FIG. 14: I'^O NMR shift calculated from Eq. ^ for U = 
1.75 for a 34 X 34 system averaged over 15 random impurity 
configurations at g^BB/2 = 0.004 and T = 0.013. Numerical 
smoothing was performed by convolving the distribution with 
a Lorentzian of width 3.5kHz. 



17 7^ / ^ _ ^''^hf [Xsptnjx +^]y)+ Xspinjx ^ ^V)] 

Ms 

(13) 

where we assumed here an oxygen lying along the x axis, 
with X and y the Cu coordinates in units of the Cu lattice 
spacing a. 

In the pure metallic state above Tc, Xspin is uniform, 
so that the NMR consists of a single line shifted by 
^'^Kspin = 2i'^A,;/Xspin/MB- In the superconducting 
state, the situation is more complex because of the pres- 
ence of the vortex lattice. Eq. H13|l stays valid, but B{y) 
is no longer uniform. As the NMR spectrum is a his- 
togram of all the frequencies i^(r), this spectrum directly 
reflects the vortex field distribution, assuming that or- 
bital and spin shifts stay uniform (see for example Ref. 

Let us now consider how the impurities affect the NMR 
spectrum. If an impurity induces a spatially dependent 
magnetization as computed above, it can be regarded as 
a distribution of the susceptibility Xspin(x]y) among Cu 
sites, leading to a distribution of Kspim i-e. a broaden- 
ing of the NMR spectrum. When both impurities and 
superconductivity are now taken into account, the NMR 
spectrum should be distributed simultaneously by a dis- 
tribution of the local field B{r) due to the vortex lattice 
and by a distribution of spin shifts Kspini'^) due to the 
impurities. However, these two distributions are uncor- 
related, as argued in Ref. 0. They should then simply 
convolve with each other. 

In Fig. 1141 we show how the distribution of magnetiza- 
tions is transformed into a distribution of ^^O frequency 
shifts due to disorder alone. As expected, increasing dis- 
order broadens the line. In addition, however, it is seen 



that the line asymmetry increases, with enhanced weight 
on the positive side. This effect was indeed observed in 
experiment, and attributed to the enhanced spin suscep- 
tibility, i.e. density of states near the Fermi level due 
to impurities. Thus some small line asymmetry would 
result simply because the density of states of a disor- 
dered d-wave superconductor is enhanced in the intersti- 
tial regions far from the impurity, but as we have seen the 
susceptibility is also enhanced due to correlation effects, 
which are also affected by interference of many impuri- 
ties. 

In order to compare our computation to the experi- 
mental NMR lines in the presence of Zn below Tc, we 
must now estimate the effects of the vortex lattice. The 
idea is to produce a series of NMR lines corresponding 
to the pure system in the presence of the vortex lattice 
field modulation, renormalized by the disorder-enhanced 
magnetic penetration depth for each impurity concentra- 
tion. These lines include in principle the effects of spa- 
tially varying superflow, but not of spatially varying spin 
magnetization. 

Simulating the vortex field distribution including 
quasiparticle contributions with high accuracy is prob- 
ably difficult even in the nominally pure case due to un- 
certainties regarding the origin and statistical nature of 
the disorder in the vortex lattice. Following Ref. 0, we 
therefore identify empirically the vortex-induced part of 
the field distribution in the superconducting state with 
the field- independent part of the overall NMR shift. We 
accept the determination of the width and shift of the 
field-independent part of the distribution determined by 
Ouazi et ai.O], shown in Fig. ^1 and estimate the vortex 
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Planar Zn concentration 



FIG. 15: Field-independent experimental ^^O shifts (squares) 
and widths (triangles) from Ref. ^ Note no field-independent 
shift was determined for the 1.5% sample. 
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FIG. 17; Top: normalized theoretical ^'^O NMR lines at 
T — 0.013 and g^BB/2 = 0.004 obtained by the procedure 
described in the text, for U = 1.75 and planar Zn concentra- 
tions of Xp = 0%, 1.5%, 3%, and 6%. Bottom: Experimental 
results on YBCO powders from Ref. i for the same concen- 
trations in a 7T field at 15K. 
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FIG. 16: Lorentzian vortex field distributions deduced from 
field-independent data of Ouazi et al 9j. Note the 1.5% curve 
was obtained by interpolating the shift of the 1.5% sample. 



contribution by the Lorentzians determined using these 
parameters. The lines thus obtained are shown in Fig. 
[TCI Note that the data for the 1.5% Zn sample was not 
good enough to extract the field- independent part; we 
have therefore simply interpolated hnearly between the 
pure and 3% sample. 

These lines are now convolved with the distributions 
of NMR shifts obtained from the impurity effects 
alone. This is justified because in fields of a few Tesla 
the inter- vortex distance is of the order of hundreds of A, 
whereas the typical inter-impurity distance for samples 
with per cent level Zn is tens of A. Thus, there can be 
no significant correlation between the positions of most 
of the Zn atoms and the vortices themselves. 

To obtain the corresponding ^^O shifts, we now use the 
values 7/27r = 5.772 MHz/T, Ahf =36kOe [sllll and 



^'^Korb — 0.02% which apply for the external field 
along the c crystallographic axis as in Ref. H Finally we 
convolve both the impurity and vortex field distributions. 
The lines thus obtained are plotted in Fig. 1171 and com- 
pared to the experimental results of Ouazi et a].j9| taken 
at 15K. The semi-quantitative variation of the width and 
asymmetry of the experimental lines with impurity con- 
centration are seen to be remarkably well reproduced by 
the theoretical results. In addition, the magnitudes of the 
shifts for different Zn concentrations are quite well repro- 
duced, with the exception of the 1.5% sample, where the 
shift of the vortex field distribution was effectively un- 
known. We note that no extensive fitting in parameter 
space was done, so it is quite striking that the magnitudes 
and dependence on concentration agree so quantitatively. 

To make further quantitative comparisons with the 
widths, which are the more experimentally reliable quan- 
tities, we plot in Fig. 1 181 the same normalized lines shown 
in Fig. El but with the shifts removed. It is seen that 
the theory tracks the increase in width as well as the 
overall lineshape extremely well. 
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FIG. 18: Normalized experimental "O NMR lines of Ref. 
at 15K and 7T with all shifts removed (solid lines) and 
normalized theoretical impurity NMR lines for U — 1.75 con- 
volved with vortex field distributions as described in text, 
shifts also removed (dashed lines): (a) pure (b) Xp = 1.5%, 
(c) Xp = 3.0%, and (d) Xp = 6.0%. 
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FIG. 19: Single impurity of strength Vimp = 100 indicated by 
circle at center of 34 x 34 system. Sites with magnetization 
values within ranges such as to determine NMR linewidth for 
U = 1.75, T = 0.01, 3MB-B/2=0.004. Sites with -0.0015 < 
m < —0.0001 ("!^L sites") are colored black, and those with 
0.0005 < m < 0.002 {"uh sites") are colored gray. 



B. T dependence of line asymmetry 

Oiiazi et observed an increase in the NMR line 

asymmetry, with a shift in weight towards the positive 
side, as Zn concentration was increased and/or as the 
temperature was lowered. They proposed that this phe- 
nomenon was associated with the formation of the reso- 
nant state around Zn observed by STM, increasing the 
LDOS at the Fermi level near the impurity, and thereby 
enhancing the spin susceptibility. Since the LDOS en- 
hancement is always positive, this effect selectively en- 
hances the broadening on the positive side of the line, 
provided it exceeds in magnitude the homogeneous mag- 
netization in the regions of the sample far from impuri- 
ties. In a d-wave superconductor, this is always true at 
sufhciently low T since the bare susceptibility of the clean 
system vanishes as ~ T at low temperatures. Note that 
this LDOS-based asymmetry enhancement is essentially 
the same effect to which the "LDOS-only" approaches 
mentioned above ascribe the entire enhancement of the 
local susceptibility measured in NMR. We have already 
argued that this is a small effect with regard to the over- 
all T dependence in the near field of the impurity; here 
we show that it can nevertheless play an important role 
in the T-dependent structure of the line, and in partic- 
ular the asymmetry of the linewidth, which arises from 
contributions further from the impurities. 

We find that the explanation given by Ouazi et ai.jOj 
is essentially correct, but is strongly enhanced both by 
interference between multiple impurities and by antifer- 
romagnetic correlations. We have seen that the suscep- 



tibility on the sites nearest the impurity is strongly in- 
creased as the temperature is lowered (Fig. |7)l, and that 
this effect is magnified by increasing U (Fig. (H)) . The en- 
hancement of the asymmetry of the lineshape due to the 
interference effect for fixed temperature was mentioned 
above, and illustrated in Fig. ^| We therefore antic- 
ipate that the low-T upturns in the asymmetry of the 
lineshape will be a characteristic of the present model as 
well. 

To see the origin of the asymmetry enhancement at 
low T, let us first examine the temperature dependence 
of the magnetization of those sites which actually de- 
termine the measured width. This is not a completely 
straightforward proposition, given that the magnetiza- 
tion patterns in the d-wave superconducting state are de- 
termined by a combination of normal Friedel oscillations, 
local magnetic correlations, d-wave pair correlations, and 
interference in the many-impurity case. Even for a single 
impurity, the first three effects combine to make it diffi- 
cult to specify, e.g. a given distance from the impurity 
which is important for determining the linewidth. It is 
clear that it is not the nearest- or next-nearest- neigh- 
bors which do so, but as shown in Fig. ^1 the set of 
sites actually contributing to the positive and negative 
half-widths - while indeed clustered around a range 10- 
15 lattice spacings from the impurity in the nodal tails 
of the wavefunctions, as noted above - form a more com- 
plicated pattern. 

The temperature dependence of these selected sites is 
now shown in Fig. I2UI where the influence of electronic 
correlations is also illustrated by comparing U = 1.75 
and U = 0. In the noninteracting case, we can see the 
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FIG. 20: Temperature dependence of magnetization on sites 
identified as contributing to determination of NMR linewidth, 
as specified in Fig. [T^ 



low-temperature upturn of the nearest-neighbor magne- 
tization, as discussed in Ref. [13 (note at T ^ 0, the 
noninteracting local susceptibility always — > since the 
impurity resonance sits at a finite energy for any generic 
potential). On the other hand, the upturn of the mag- 
netization of the sites contributing to the linewidth is 
much weaker to nonexistent in the noninteracting case, 
although the LDOS effect still manifests itself via the 
fact that A;/i < /S.vh- In the [/ > case, however, 
the upturns are much stronger, and even manifest the 
saturation of the magnetization arising from Aj/^ sites 
observed in experiment 0] , which leads to the surprising 
increase of the /S.vh — observed. 

To verify if the phenomenon is reproduced quantita- 
tively, we plot explicitly the difference of the half width 
at half-maximum on the high frequency side of the O 
line, /S.VH, and the same quantity on the lower side /S.vl. 
Each is independently enhanced at low temperatures, but 
it is the difference which is particularly striking, as seen 
in Fig. 1^ where theory for U = 1.75 is compared with 
the Ouazi et ai.Q results shown in Fig. [221 While the de- 
tails of the theoretical curves do not agree exactly with 
experiment, it is clear that the basic results are repro- 
duced by the theory, both in terms of the temperatures 
at which the upturns begin, and in terms of the magni- 
tudes of the upturns themselves. On the other hand, the 
lower halfwidth l^v^, is roughly T independent in exper- 
iment, but has a weaker but still significant enhancement 
in the calculation, as indicated in Fig. ^| We do not un- 
derstand the origin of this discrepancy at present. In ad- 
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FIG. 21: Line asymmetry-theory: difference of halfwidths on 
negative {^vl) and positive {/S.vh) frequency sides of NMR 
lines vs. T for five impurity concentrations. Top panel: U = 
1.75, pAtsB/2=0.004, Vi^p=im. The result for U=Q is also 
shown for comparison (dotted line). Bottom panel: same, but 
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dition, the theoretical result retains a certain asymmetry 
of the lineshape up to higher temperatures, whereas the 
experimental lineshape becomes symmetric above about 
40K. This may be due to the neglect of inelastic scatter- 
ing, which becomes important at higher T, in the calcu- 
lation. 

The two figures shown in Fig. 1211 correspond to two 
different values of the impurity potential Vimp- The up- 
per panel corresponds to a value Vimp = 100, which gives 
an impurity resonance energy of fio ~ — O.Oli, whereas 
the Vimp = 10 lower panel corresponds to a resonance 
energy of JIq — within our numerical resolution, for 
this particular band. Thus the local Fermi level density 
of states near the impurity is larger in the second case, 
leading indeed to a stronger upturn, as anticipated. The 
actual resonance energy of a Zn in YBCO is not known 
at this time, but is expected to be close to the -l.SmeV 
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FIG. 22; Line asymmetry-experiment: difference of 
halfwidths on negative {^vl) and positive {Avh) frequency 
sides of NMR lines vs. T for three impurity concentrations, 
from Ref. H 
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50 



T(K) 



100 



150 



200 



Q. 



1000 



500 











1.0% 






-- 2.0% 






3.0% 






- ■- 6.0% 















0.05 



0.1 
T 



0.15 



0.2 




FIG. 23: Upper panel: Li Knight shift as calculated in text, 
using same magnetizations, obtained for same parameters as 
shown in the histograms in Fig. 1111 3 random configurations. 
Lower panel: experimental Knight shift data from Ref. ^ for 
Xp^ 1% and 2% Li in optimally doped YBCO for B=7T. 



V. FURTHER COMPARISONS WITH 
EXPERIMENT: ^LI 

The "universality" of the magnetic response to strong 
in-plane nonmagnetic defects was noted early on by Bo- 
broff et al. Q]. That is, any in-plane impurity which 
from a conventional chemistry standpoint is expected to 
be nonmagnetic appears to have a nearly identical effect 
on both normal state and superconducting properties. 
This includes Zn, Li and defects in the plane created 
by electron irradiation, which produce nearly identical 
changes in susceptibility, Tc, and resistivity per in-plane 
impurity. This is remarkable because there would appear 
to be important electronic structure differences between 
the Zn ion, which has a closed shell, and a Li, which is 
believed to localize a hole around itself. It is generally 
believed that the essential features of these defects are 
therefore simply their ability to exclude mobile conduc- 
tion electrons, hence our choice of model of, e.g. Zn as 
a strong repulsive potential. Within this assumption of 
"universality" of in-plane nonmagnetic defects, we can 
take the results for the magnetization distribution in the 
disordered system already produced, and use them to 
describe the results of earlier ^Li NMR experiments on 
YBCO. 

Li has the thus far unique ability in the cuprates to 
simultaneously provide an in-plane impurity and a nu- 
cleus (^Li) suitable for NMR. The signal is therefore not 
complicated by contributions from regions of the sample 
far from impurity sites, but provides direct information 
about the immediate vicinity of the impurity, which re- 
places a Cu in the Cu02 plane. We will assume, as in 



prior work, that the Li provides a shift equal to the sum 
of the magnetizations on its nearest neighbor sites, lead- 
ing to a shift (compare Eq. 1)13(1 '). 

^g.,.„(r)^ '^^^^^^-^-^" + ^\ (14) 

fJ-B 

where 5 is a nearest neighbor displacement. Results plot- 
ted as a function of temperature are shown in Fig. 1 2 31 and 
compared to the experimental results on Li substituted in 
optimally doped YBCO of Bobroff et al. 8]. The theory 
seems to reproduce the initial weak dependence on the 
Li concentration (nonexistent in the normal state, small 
in the superconducting state). This is because the Li are 
sensitive to the local environment of the impurity only, 
and the shift does not therefore depend on concentration 
until interference effects become significant. Until now, 
experiments have only been performed on a maximum 
of 2% planar Li impurities; the deviations of the curves 
corresponding to higher concentrations at low tempera- 
tures constitute predictions of the theory which can be 
verified by NMR. We note that the experimental results, 
in particular the change in slope near Tc and the sub- 
sequent rise at low T, have been discussed in terms of 
the effect of the opening of the superconducting gap on 
Kondo screening of simple moments by the d-wave quasi- 
particle gas|^. Our results here, which do not account 
for the spin-flip scattering necessary for the Kondo effect 
to take place at all, indicate that a strong-correlation ex- 
planation for these, and associated phenomena, is more 
likely. 
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VI. CONCLUSIONS 

We have shown that a theoretical model of strong im- 
purities in a d-wave superconductor, with residual quasi- 
particle interactions treated in a weak coupling Hub- 
bard model within a mean-field approximation, provides 
an excellent description of NMR experiments on opti- 
mally doped YBCO. This model assumes that both Zn 
or Li impurities act as simple strong potential scatter- 
ers which induce staggered magnetic polarization clouds 
around themselves with very small total moment pro- 
portional to the external magnetic field. The results and 
comparison with experiments on optimally doped YBCO 
are not compatible with static impurity-induced mag- 
netism present in zero field. This is consistent with the 
lack of a magnetic /iSR signal over most of the YBCO 
phase diagram down to very low doping [Tsi? ], and 
with the absence of a Schottky anomaly in specific heat 
measurements |i38l| . The universal response of the system 
to nonmagnetic defects which are very different chemi- 
cally also confirms the primacy of the strong correlations 
in the Cu02 plane for these phenomena. 

The ability of NMR to probe nuclei with different con- 
figurations with respect to the Cu spins means that dif- 
ferent interesting aspects of the physics can be tested by 
the model. In particular, the ^^O line is primarily a probe 
of the weaker magnetism induced far from the impurities. 
It is thus more sensitive to the effects of interference of 
the quasiparticle bound states associated with different 
scattering centers. We have shown that this quasiparticle 
interference is responsible for the enhancement of the ^^O 
linewidth as measured by experiment, and the variation 
of this lineshape with Zn concentration has been repro- 
duced quantitatively by the model. The low-temperature 
enhancement of the ^^O line to the positive side caused 
by Zn appears, within this framework, to be related to 
the LDOS enhancement around the impurities, as sug- 
gested by Ouazi et aI0- On the other hand, the size 
and strong temperature dependence of the line broaden- 
ing can be understood only if electronic correlations and 
concomitant paramagnetic local moments are present. 

The ^Li nucleus, which senses only the nearest neigh- 
bor Cu, gives us a picture into magnetic effects in the 
vicinity of the individual impurities. Our study, which 
reproduces quantitatively the ""Li Knight shift magnitude 
and T dependence in both the normal and superconduct- 
ing states, concludes that experiments on Li done thus 
far have probed only the regime of individual isolated im- 
purities, where interference effects are not strong enough 
to affect nearest neighbor sites, but these should be visi- 
ble by going to only slightly higher Li concentrations. In 
addition, it shows how resonant states in the d-wave su- 
perconductor enhance both the single-impurity magnetic 
effects and the quasiparticle interference. 

We note a number of questions which are open and 
which should be clarified in subsequent work. First, 



the current framework is a weak-coupling mean field ap- 
proach which neglects the renormalization of the elec- 
tronic structure near the impurity sites found in strong- 
coupling approaches. We anticipate that such effects, 
which arise from diagrams not included in the RPA-type 
analysis here, will give rise to quantitative changes in the 
values of parameters chosen here to fit experiment, but 
not alter the overall physical picture we have presented. 
Many of these issues have already been raised and un- 
derstood in the context of homogeneous RPA-level spin- 
fluctuation theories, where a reduced effective interaction 
U replaces the bare U after resummation of a subclass of 
diagrams. It is intriguing to note that NMR experiments 
seem to require us to work very close to the magnetic 
phase transition in the theory in order to explain e.g. 
the magnitude of the linewidth and the T dependence of 
the line asymmetries. This is reminiscent of the early fits 
of spin fluctuation theories of NMR for the homogeneous 
systems, which required values of the effective interaction 
to be close to the bulk transition in the model. Taken 
together, it is tempting to speculate that these results 
indicate the flow of the true Hamiltonian of the system 
to strong coupling at sufficiently low energies. 

In this work we have not investigated carefully the 
transition to the normal state, which we reserve for a 
subsequent paper. It is interesting to note that, while ear- 
lier works on the normal state using the same model0,0 
found it necessary to introduce an extended impurity po- 
tential to fit experiment, our conclusion is that this is 
not necessary in the superconducting state. Following 
Bulut's suggestion^J that the extended nature of scat- 
tering potential may arise in the normal state from the 
renormalization of the local electronic structure by cor- 
relations, this may be an indication that the processes 
leading to this effect are simply gapped below Tc- 

Finally, we emphasize that our calculation does not 
include spin-flip scattering terms necessary to recover a 
true Kondo effect, and that to the extent our analysis has 
been successful, our results therefore imply that Kondo 
physics is probably unnecessary to describe the phenom- 
ena in question, at least at optimal doping. It will be 
interesting to see whether the phenomenology put for- 
ward here continues to hold as one goes to underdoped 
systems with planar impurities which remain paramag- 
netic. Other possible extensions include more strongly 
interacting and/or more disordered systems, where static 
magnetism in zero field may be present, studies of the 
normal state, as well as of the overdoped regime. Work 
along these lines is in progress. 
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